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摘要 : 园 环 结构 的 准确 识别 和 跟踪 有 利于 更 好 地 研究 和 分 析 园 环 特 征 以 及 同 光斑 爆发 和 
日 园 物 质 抛射 之 间 的 关系 。 提 出 基于 相位 一 致 性 和 方向 滤 疲 的 园 环 震荡 识别 和 跟踪 方法 。 识 
别 过 程 如 下 首先， 使 用 相位 一 致 性 技术 对 太阳 图 像 中 的 园 环 特征 进行 识别 ; 其 次 ， 对 识别 
后 的 图 像 进 行 二 值 化 处 理 ， 接 着 使 用 方向 滤波 去 除非 园 环 的 结构 特征 ;， 然后， 对 其 进行 形态 
学 处 理 ; 最 后 ， 再 对 识别 结果 进行 二 次 曲线 拟 合 。 运 用 所 提 算 法 对 2011 年 9 月 6 日 由 太阳 
动力 学 观测 站 (Solar Dynamics Observatory, SDO) Æ Fe IX 17. 1 nm 波段 观测 到 的 园 环 震荡 进 
行 了 识别 和 跟踪 ， 结 果 验 证 了 所 提 算 法 的 可 行 性 和 准确 性 。 证 明 所 提 算 法 能 被 用 于 冕 环 特 性 
以 及 震荡 过 程 的 研究 和 分 析 。 
关键 词 : 相位 一 致 性 ; 方向 滤波 ; ARED; 震荡 周期 
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日 园 位 于 太阳 大 气 的 最 外 层 ， 在 上 面 可 以 观测 到 一 些 亮 的 环形 结构 ， 它 们 被 称 为 园 环 。 园 环 由 日 
冕 磁场 束缚 的 热 等 离子 体 发 亮 而 形成 环 状 结构 ， 是 太阳 磁场 结构 的 一 种 表现 形式 。 它 们 同日 冕 物质 抛 
射 、 太 阳 耀 班 、 太 阳 风 暴 等 有 密切 的 关系 。 研 究 冕 环 及 其 演化 过 程 对 于 理解 太阳 大 气 层 的 磁场 活 
动 和 预报 太阳 风暴 具有 重要 意义 。 

太阳 溜 斑 爆发 和 日 冕 物质 抛射 释放 的 大 量 能 量 引起 附近 局 部 磁场 结构 的 急剧 变化 ， 导 致 附近 的 岗 
环 结构 产生 震荡 。 这 些 震 荡 不 仅仅 发 生 在 一 条 冕 环 结构 上 ， 也 同时 出 现在 多 条 冕 环 结构 上 。 最 近 
的 研究 表明 ， 这 种 震荡 是 驻 波 的 一 种 表现 形式 ， 但 是 每 个 环 上 由 于 震荡 的 频率 、 相 位 以 及 衰减 系数 的 
不 同 ， 往 往 呈 现 不 同 的 震荡 特征 。 研 究 这 些 特征 对 于 更 好 地 理解 炊 斑 爆发 和 日 园 物 质 抛 射 对 日 园 
人 磁场 的 影响 具有 很 大 的 帮助 。 

然而 准确 地 识别 和 提取 冕 环 特征 结构 面临 非常 大 的 困难 。 这 是 由 于 : (1) 冕 环 是 由 日 冕 磁场 束缚 
的 热 等 离子 体形 成 ， 这 些 等 离子 体 分 布 并 不 均匀 ， 虽 然 晃 环 结构 的 强度 较 大 ， 但 强度 变化 (梯度 ) 较 
小 ， 导 致 野 环 结构 的 边缘 模糊 ， 造 成 提取 困难 ; 〈2) 各 个 网 环 之 间 的 距离 较 近 而 且 相 互 交 错 ， 导 致 产 
生 的 二 维 图 像 环 形 结构 彼此 交叉 重大 ; (3) 由 于 磁场 活动 的 复杂 性 导致 了 网 环 结构 呈现 多 种 形态 ， 如 
环形 结构 出 现 断 裂 等 问题 ; (4) 在 冕 环 附近 还 存在 非 环 形 的 太阳 特征 ， 如 海绵 状 亮 班 ， 这 些 都 增加 了 
冕 环 提取 的 复杂 性 和 准确 性 。 

目前 冕 环 的 识别 方法 主要 分 为 两 类 "“” : 一 类 是 运用 图 像 空间 域 的 强度 和 梯度 阔 值 对 特征 结构 
进行 识别 和 提取 ; 男 一 类 是 利用 频 域 或 小 波 域 对 不 同 的 频率 特征 进行 滤波 ,识别 和 提取 园 环 特征 。 
广义 上 讲 ， 图 像 的 结构 特征 ， 如 特征 边缘 、 尺 度 和 纹理 等 大 量 底层 信息 存储 在 图 像 的 相位 谱 中 。 文 
[14-15] 证 明了 相位 一 致 性 在 表征 特征 结构 时 对 强度 、 对 比 度 、 照 度 具 有 不 变性 ， 而 且 通 过 计算 图 像 
的 局 部 能 量 实 现 了 对 图 像 特征 的 稳定 提取 。 目 前 ， 相 位 一 任性 技术 已 经 应 用 于 黑子 本 影 亮点 和 本 影 炊 
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斑 以 及 光 球 米粒 等 低 对 比 度 太 阳 图 像 特征 的 提取 和 识别 “” 。 
本 文 运用 相位 一 致 性 技术 并 结合 方向 滤波 的 方法 对 一 个 由 太阳 动力 学 观测 站 搭载 的 太阳 大 气 成 像 
仪 (Atmospheric Imaging Assembly, AIA ) Æ 2011 年 9 月 6 日 由 谱 线 Fe IX 17. 1 nm PEER WLI BY AH 
震荡 进行 识别 和 分 析 ， 并 同 已 有 的 文献 结果 进行 比较 来 检验 所 提 方 法 的 可 行 性 和 准确 性 。 
论文 首先 介绍 所 用 的 实验 数据 ， 其 次 们 述 早 环 震荡 的 识别 过 程 ， 然 后 对 识别 的 蝎 环 震荡 进行 统计 
分 析 ， 并 同 已 有 文献 结果 作 比 较 ， 最 后 对 本 文 进行 总 结 。 


1 实验 数据 


选择 震荡 事件 发 生 在 编号 为 11283 的 活动 区 为 实验 数据 ， 时 间 从 2011 年 9 月 6 日 22 时 17 分 到 2011 
年 9 月 6 日 22 时 34 分。 观测 资料 由 太阳 大 气 成 像 仪 在 谱 线 Fe IX 17. 1 nm 波段 记录 。 园 环 震荡 由 一 个 附 
近 爆 发 的 X2. 1 级 次 斑 引 起 。 该 组 数据 共 包 含 65 幅 图 像 ， 间 隔 时 间 为 12 s， 图 像 的 像 元 分 辩 率 是 0.6", 
图 1(a) 展 示 了 2011 年 9 月 6 日 22 时 17 分 获得 的 图 像 其 中 白色 稍 头 所 示 的 位 置 为 次 斑 爆发 的 位 置 ， 
它 引 起 了 图 1(a) 中 白色 方 框 标记 区 ( 即 图 1(b) ) 园 环 产 生 了 震荡 。 图 1(b) 展 示 了 这 个 区 域 的 细节 。 它 的 
大 小 为 210 x 144"。 在 图 1(b) 中 ， 由 数字 1 和 2 指示 的 位 置 为 次 斑 爆 发 能 量 释放 所 产生 的 图 像 结构 特征 ， 
它 导 致 了 附近 的 园 环 产生 剧烈 震荡 ， 尤 其 是 第 头 Loop 4 和 Loop B 指示 的 两 个 冕 环 结构 震荡 尤为 明显 。 
为 了 同文 [9] 做 比较 ， 本 文 主要 针对 冕 环 4 AID B 的 震荡 运用 所 提 算 法 进行 识别 和 跟踪 ， 并 对 其 结 
果 进 行 统计 分 析 来 检验 算法 的 可 行 性 和 准确 性 。 
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图 1 ZA 
Fig. 1 Loops 


2 识别 和 跟踪 方法 


算法 首先 运用 文 [14-15] 提 出 的 图 像 的 相位 一 致 性 方法 计算 图 1(b) 中 的 相位 一 致 性 特征 。 图 2 
(a) 是 获得 相位 一 致 性 中 的 局 部 加 权 平 均 相 位 角 图 像 。 它 的 取 值 范围 在 [-m/2, 5/2] 之 间 。-m/2 
和 TZ2 分 别 表示 原始 图 像 中 一 个 局 部 区 域内 明暗 特征 极 值 的 位 置 。0 表示 特征 由 明 到 暗 (或 由 暗 到 
明 ) 的 分 界 。 从 图 2(a) 可 以 看 到 , 在 图 1(b) 中 的 环形 结构 以 及 由 炊 斑 爆发 产生 的 特征 和 海绵 状 亮 斑 
也 被 同时 识别 出 来 。 由 于 在 原始 图 像 中 冕 环 结构 为 亮度 较 大 的 特征 ( 即 冕 环比 较 明 亮 ) ， 因 此 ， 识 别 
算法 在 对 图 2( a) 进行 二 值 化 处 理 时 选择 阀 值 为 0， 即 图 2(a) 中 大 于 0 的 被 置 为 1， 和 否则 为 0， 其 二 值 
化 结果 如 图 2(b)。 

从 图 2(b) 可 以 看 出 ， 那 些 海绵 状 亮 斑 和 由 于 光斑 爆发 产生 的 近似 垂直 的 一 些 特征 (如 图 中 上 
方 中 间 位 置 被 白色 方 框 标记 ) 干 扰 冕 环 4 和 下 的 提取 。 所 要 识别 和 提取 的 冕 环 4 ASIA BOT HE 
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似 水 平 。 因 此 算法 选用 方向 滤波 思想 识别 近似 垂直 的 特征 并 去 除 。 算 法 选用 sobel 边缘 检测 算 子 
g(g=[1, 0, -1; 2, 0, -2; 1, 0, -1]) 作 为 方向 滤波 算 子 检测 干扰 线段 的 垂直 边缘 并 去 除 它 们 ， 处 理 
结果 如 图 2(c)。 

从 图 2(c) 可 以 发 现 , 那些 近似 垂直 的 干扰 特征 通过 使 用 方向 滤波 的 方式 已 经 被 很 准确 地 去 除 。 
然而 所 需 提取 的 网 环 结构 仍然 存在 不 连续 以 及 冕 环 上 存在 小 的 孔洞 且 边 缘 不 规则 、 噪 声 点 大 和 网 环 粗 
细 不 均匀 等 问题 。 算 法 采用 形态 学 的 去 除 小 面积 、 填 充 、 脱 胀 、 闭 运算 、 细 化 等 一 系列 操作 对 上 述 问 
题 进行 处 理 ， 处 理 结果 如 图 2(d) 。 

从 图 2(d) 箭 头 处 可 以 看 出 ， 识 别 的 冕 环 结构 不 光滑 ， 存 在 许多 分 支 结 构 ， 并 且 图 中 除了 所 要 识 
HNN BH A 和 冕 环 有 外 ， 还 有 一 些 短 的 非 目 标 线 。 因 此 算法 根据 节点 处 各 分 校 的 长 短 ， 保 留 最 长 的 
分 村 结构 ， 去 除 短 的 分 校 杂 点 ， 并 根据 环 的 长 度 ， 提 取 了 园 环 4 BM B, ZARA 2(e). 

最 后 对 所 识别 的 65 AVR PHBH A AH B 进行 了 二 次 方程 的 曲线 拟 合 ， 获 取 一 个 光滑 的 弧 
线 ， 其 二 次 方程 定义 如 下 : 


y=ax thx te, (1) 


其 中 ,x WERP EARME; y 为 图 像 中 冕 环 的 纵 坐 标 ; a. b, e 为 参数 。 图 2(e) PRAHA ME 
环 B 二 次 拟 合 参 数 a、b、c 的 值 如 表 1。 图 2(e) 环 的 拟 合 结果 被 芭 加 到 原 图 上 并 显示 在 图 2(f) 中。 
Hh Fee a AP EA A FH B 的 形态 均 不 相同 ， 因 此 本 文 分 别 对 其 进行 了 逐一 拟 
合 ， 其 它 的 拟 合 参数 并 未 在 此 列 出 ， 这 些 拟 合 线 为 本 文 要 进行 统计 分 析 的 特征 。 


(e) 


图 2 识别 过 程 
Fig.2 The process of Identification 


3 实验 结果 与 分 析 #1 对 图 1(b) HERA 和 用 的 二 次 拟 合 参数 
Table 1 Parameters of quadratic fitting 

WT SOE EN BA, SPER a ; p 
4 REI B EE X FTE 10 个 位 置 ， 然 后 根 BRA 0.002801. L 088 19.33 
据 每 幅 图 的 二 次 拟 合 方程 (1) 得 到 了 值 ， 并 对 这 nee ee ee 


10 个 了 值 分 别 取 它们 在 每 1 幅 图 像 上 的 均值 作 

为 这 一 幅 图 像 晃 环 震 荡 的 位 置 。 图 3(a) 和 (b) 中 星 号 标记 的 位 置 即 为 冕 环 4 和 园 环 B 的 了 均值 位 
置 随时 间 了 的 变化 情况 。 在 图 3 中 ， 横 轴 表 示 时 间 ( 单 位 为 分 钟 ) ， 纵 轴 表 示 冕 环 移 动 距离 (单位 为 
千 米 )。 从 图 中 可 以 看 出 这 些 星 号 呈现 一 个 类 似 衰减 的 正弦 函数 ， 因 此 采用 如 下 子 数 拟 合 这 些 位 移 ， 具 
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体 的 衰减 正弦 拟 合 函 数 定义 如 下 : 


f(t) = 4e "sin(27t/T + 7e/180), (2) 
其 中 ,4 为 振幅 ; 7 为 衰减 时 间 ; 上 为 时 间 ( 单 位 为 分 钟 ) ; 了 为 周期 ; gp 是 初 相位 (单位 为 度 ) 。 
1 000 
800 * Distance vs. Time * Distance vs. Time 
Loop A Oscillation — Loop B Oscillation 
600 
400 
Ss op S 
A -200 A 
-400 
-600 4 
-800 
Time/min Time/min 
(a) (b) 
图 3 震荡 
Fig.3 Oscillations 
= 对 于 这 两 个 曲线 的 震荡 周期 ， 是 采用 傅 里 叶 变换 的 频率 分 析 方 法 获得 。 这 两 个 振荡 信和 号 的 频谱 分 
一 析 图 如 图 4。 依据 伟 里 叶 频谱 分 析 ， 这 两 个 环 的 震荡 频率 约 为 8.4 MHz， 周 期 为 118 s。 运 用 这 个 周期 


值 T( HU 118 s) 和 拟 合 方程 (2) 拟 合 图 3 中 的 星 号 位 置 。 拟 合 的 曲线 如 图 3(a) 和 (D) ， 它 们 的 衰减 系数 
约 为 122s， 但 两 个 震 功 序 列 的 相位 差 近似 41" ， 冕 环 刀 先 于 冕 环 4， 而 且 有 环 的 振幅 也 高 于 4 环 。 这 
N ZO ETT RAGA EAE LEM BIE, TES BBL EN BES. 
> 不 过 这 种 位 置 的 远近 也 可 能 是 由 成 像 的 视角 导致 的 。 这 个 分 析 结 果 同文 [9] 的 结论 一 致 ， 进 一 步 验证 
了 本 文 所 提 的 方法 可 以 运用 于 网 环 特征 的 研究 和 分 析 。 


Amplitude value 
Amplitude value 


0 0L_ ae 
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 
Frequency/Hz Frequency/Hz 
(a) (b) 


图 4 频谱 分 析 图 


Fig.4 Spectrum analysis diagram 


4 讨论 与 总 结 


大 多 数 的 网 环 震 沪 是 由 附近 的 兆 斑 爆发 和 日 哆 物质 抛射 引起 的 ， 对 这 些 冕 环 的 震荡 进行 分 析 可 以 
更 好 地 理解 冕 环 的 磁场 特征 与 太阳 耀 斑 爆发 和 日 冕 物质 抛射 的 能 量 释放 之 间 的 关系 ， 准 确 地 识别 和 提 
取 冕 环 特征 至 关 重 要 。 

本 文 运用 相位 一 致 性 和 方向 滤波 的 方法 分 析 了 一 个 由 太阳 大 气 成 像 仪 在 17. 1 nm 波段 观测 到 的 冕 
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RH. AT, TRIIN, KET Ee KRETA HE RL 

致 性 和 方向 滤波 的 图 像 特征 提取 方法 识别 和 追踪 它们 ， 并 根据 识别 和 追踪 的 结果 对 这 两 条 冕 环 的 周期 、 

相位 差 以 及 衰减 系数 进行 统计 分 析 ， 并 同 已 有 文献 的 结果 进行 比较 ， 证 明 算 法 可 行 ， 识 别 结果 准确 。 
尽管 相位 一 致 性 可 以 很 好 地 解决 低 信 噪 比 条 件 下 图 像 特征 提取 的 问题 ， 但 由 于 网 环 特征 存在 于 复 


dab E. 
杂 的 背景 


之 中 以 及 冕 环 结构 边缘 模糊 和 信 品 比 低 等 因素 ， 使 得 提取 的 冕 环 结构 不 太 完整 ， 一 些 比较 


细 、 瞳 的 晃 环 结构 未 能 识别 成 功 。 对 于 这 些 情况 仅仅 从 图 像 处理 的 角度 识别 和 提取 是 不 够 的 ， 需 要 结 
合 其 他 波段 的 数据 及 磁场 反 演 等 物理 过 程 分 析 ， 提 高 蝎 环 识别 的 准确 性 ， 这 也 是 下 一 阶段 研究 工作 的 


重点 。 
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Abstract; Solar flares occurring usually cause oscillations of nearby coronal loops. A precondition of 
studying these oscillations is to explore how to identify coronal loops. Meanwhile, accurate identification is also 
crucial to study and analyze the statistical features of the coronal loops and the relationship to coronal magnetic 
field. However, the traditional identification methods based on intensity and/or gradient threshold have some 
limitations, because coronal loops hold a relatively low contrast, blurred boundary, and the image noise is 
close to coronal loops. So we propose a new method to identify coronal loops using phase congruency and 
directional filter. The identification procedures are as follows; firstly, phase congruency technique is used to 
identify linear shape in a coronal image; secondly, a binary image is obtained from the image of phase 
congruency for detecting the loop features; thirdly, a directional filter based on Sobel Operator is applied to 
remove the redundant features that are nearly perpendicular to those target loops; fourthly, morphological 
filters are used to extract the shape of those loops; finally, quadratic fitting is used to smooth the identified 
loop features. For illustrating the performance of the method, the coronal loop oscillations taken with the 
Atmospheric Imaging Assembly on aboard the Solar Dynamic Observatory on September 11, 2011 are used to 
evaluate the identification and tracking process. The coronal loop oscillations were caused by a nearby X2. 1 
flare emission. We analyzed the identified results, and calculated their oscillation period, phase difference and 
decay coefficient. They are in good agreement with those of previous works. This implies that the proposed 
method can accurately identify the shape of a loop in an image; furthermore, the results can also be used to 
study the physical properties of coronal loop oscillations and the relationship to others solar phenomena. 
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